arXiv: 1501.0485Ivl [gr-qc] 20 Jan 2015 


Introduction to the application of the dynamical systems theory in the study of the 
dynamics of cosmological models of dark energy 

Ricardo Garcia-Salcedo,^’EI Tame Gonzalez,Francisco A. 

Horta-Rangel, ^’El Israel Quiros,^’EI and Daniel Sancliez-Guzman^’0 

^ CICATA - Legaria del Instituto Politecnico Nacional, 11500, Mexico, D.F., Mexico. 

^Departamento de Ingenieria Civil, Division de Ingenieria, Universidad de Guanajuato, Gto., Mexico. 

(Dated: January 21, 2015) 

The theory of the dynamical systems is a very complex subject which has brought several surprises 
in the recent past in connection with the theory of chaos and fractals. The application of the 
tools of the dynamical systems in cosmological settings is less known in spite of the amount of 
published scientific papers on this subject. In this paper a - mostly pedagogical - introduction 
to the application in cosmology of the basic tools of the dynamical systems theory is presented. 

It is shown that, in spite of their amazing simplicity, these allow to extract essential information 
on the asymptotic dynamics of a wide variety of cosmological models. The power of these tools 
is illustrated within the context of the so called ACDM and scalar field models of dark energy. 

This paper is suitable for teachers, undergraduate and postgraduate students from physics and 
mathematics disciplines. 
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I. INTRODUCTION 


Gosmology is about the study of the structure and evo¬ 
lution of the universe as a whole. Einstein’s general rel¬ 
ativity (GR) stands as its mathematical basis. In gen¬ 
eral, Einstein’s GR equations are a very complex sys¬ 
tem of coupled non linear differential equations, but it is 
simplified by the underlying spacetime symmetries. The 
Friedmann-Robertson-Walker (FRW) spacetime model is 
depicted by the line element^ 


= —dt^ + a{ty 


dr^ 


1 — kr'-‘ 


-b -ksin^ Odcj)^) 


where t is the cosmic time, r, 6 and (j>, are co-moving 
spherical coordinates, a{t) is the cosmological scale fac¬ 
tor, and the constant fc = —1,0, -|-1 parametrizes the cur¬ 
vature of the spatial sections, is based on what is known 
as the cosmological principle: in the large scales the Uni¬ 
verse is homogeneous and isotropic. Homogeneity and 
isotropy imply that there is no preferred place and no 
preferred direction in the universe.^ The observations 
confirm the validity of the cosmological principle at large 
scales 100 Mpc ~ 10^® cm). The dynamics of the uni¬ 
verse is fully given by the explicit form of the scale factor 
which depends on the symmetry properties and on the 
matter content of the Universe. 
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^ For a pedagogical exposition of the geometrical and physical sig¬ 
nificance of the FRW metric see Q. 

^ Here we adopt the units where 8itG = c = 1. 


In order to determine the scale factor it is necessary to 
solve the cosmological Einstein’s field equations, which 
relate the spacetime geometry with the distribution of 
matter in the Universe. The maximal symmetry implied 
by the cosmological principle allows for great simplifica¬ 
tion of the Einstein’s equations. In this paper, for further 
simplification, we shall consider FRW spacetimes with 
flat spatial sections (fc = 0). The resulting cosmological 
equations are written as it follows:^ 


31?^ ='^Px + Px, 

X 

^ = ~\ i^lxPx+lXP^ , 

Px-kilxHpx = ^, Px+^lxHpx (1) 

where H = d/a is the Hubble parameter, px and px are 
the energy densities of the gravitating matter and of the 
A-fluid (the dark energy in the present paper) respec¬ 
tively, and the sum is over all of the gravitating mat¬ 
ter species living in the FRW spacetime (dark matter, 
baryons, radiation, etc.). The last two equations above 
are the conservation equations for the gravitating matter 
degrees of freedom and for the A-fluid respectively."^ Be¬ 
sides, it is adopted that the following equations of state 
are obeyed: 


Px = {lx - ^)Px, Px = {lx - l)pjf, 


^ For a didactic derivation of the cosmological equations not based 
on general relativity see 

^ It should be pointed out that one equation in 0, say the second 
(Raychaudhuri) equation, is redundant. 
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where Px is the barotropic pressure of the x-matter fluid, 
Px is the parametric pressure of X-component,^ while 7 ^, 
and 7 x are their barotropic parameters.® 

According to our current understanding of the expan¬ 
sion history, the Universe was born out of a set of initial 
conditions known as “big-bang.” Further evolution led to 
a stage of primeval inflation driven by a scalar field 
The inflation era was followed by a stage dominated by 
the radiation, followed, in turn, by an intermediate stage 
dominated by non-relativistic cold dark matter (CDM).^ 
During this stage of the cosmic expansion most of the 
structure we observe was formed. The end of the matter- 
domination era can be traced back to a recent moment 
of the cosmic expansion where a very peculiar form of 
matter which does not interact with baryons, radiation 
or any other form of “visible” matter - known as “dark 
energy” (DE) - started dominating Q. This alien form 
of matter antigravitates and causes the present Universe 
to expand at an accelerating pace instead of decelerating, 
as it would for a matter - either CDM or baryons - or ra¬ 
diation dominated Universe. Actually, if one substitutes 


a 

back into the second equation in dm, and taking into 
account the Friedmann equation - first equation in CD- 
then, one is left with the following equation; 


It is known that for standard (gravitating) forms of 
matter the barotropic pressure is a non-negative quan¬ 
tity: Px = hx - ^)px > 0 ^ 7a; > 1- In particular, 
for CDM and baryons the (constant) barotropic index 
7 m = 1, while for radiation jr = 4/3. Hence, since 
obviously 7 ^, > 2/3 - if forget for a while about the 
second term in the RHS of ([2) - it is seen that a cos¬ 
mic background composed of standard matter will ex¬ 
pand at a decelerated pace (d < 0). In contrast, for 
unconventional matter with ■jx < 2/3 (the parametric 
pressure px is obviously negative), the second term in 
the RHS of Eq. ([2 contributes towards accelerating 


® By parametric pressure we understand that it might not corre- 
spond to standard barotropic pressure with the usual thermo¬ 
dynamic properties but that it is rather a conveniently defined 
parameter which obeys similar equations than its matter coun¬ 
terpart. 

® In several places in this paper we use the so called equation of 
state (EOS) parameter lox = PxIpx = 7x ~ 1? interchangeably 
with the barotropic parameter 7x • 

^ Dark matter is a conventional form of matter in the sense that 
it gravitates as any other known sort of standard matter such 
as radiation, baryons, etc. However, unlike the latter forms of 
matter, CDM does not interact with radiation. This is why it is 
called as “dark matter.” 


the pace of the cosmological expansion. In particular, 
if we choose the X-component in the form of vacuum en- 
ergy (7vac,Pvac) {lx,Px), usually identified with the 
energy density of the cosmological constant (see section 
UnD, since Pvac = -Pvac, then 7 vac = 0 < 2/3. Due to 
the known fact that the vacuum energy is not diluted by 
the cosmic expansion,® eventually this component of the 
cosmic budget will dominate the cosmological dynamics, 
leading to late-time acceleration of the expansion 


a d 

A related (dimensionless) parameter which measures 
whether the expansion is accelerated or decelerated, is 
the so called deceleration parameter which is defined as 


q = 




(3) 


The pace of the expansion accelerates if g < 0. On the 
contrary, positive q > 0 means that the expansion is de¬ 
celerating. This parameter will be useful in the discussion 
in the next sections. 

In spite of the apparent simplicity of the equations CD; 
as with any system of non-linear second-order differential 
equations, it is a very difficult (and perhaps unsuccess¬ 
ful) task to find exact solutions. Even when an analytic 
solution can be found it will not be unique but just one in 
a large set of them @ . This is in addition to the question 
about stability of given solutions. In this case the dynam¬ 
ical systems tools come to our rescue. These very simple 
tools give us the possibility to correlate such important 
concepts like past and future attractors (also saddle equi¬ 
librium points) in the phase space, with generic solutions 
of the set of equations CD without the need to analyti¬ 
cally solve them. 

In correspondence with the above mentioned stages of 
the cosmic expansion, one expects that the state space 
(also, phase space) of any feasible cosmological model 
should be characterized by - at least one of - the fol¬ 
lowing equilibrium configurations in the space of cos¬ 
mological states: (i) a scalar field dominated (inflation¬ 
ary) past attractor, (ii) saddle critical points associated 
with radiation-domination, radiation-matter scaling, and 
matter-domination, and either (iii) a matter/DE scaling 
or (iv) a de Sitter® future attractors. Any heteroclinic 


® We recall that, in contrast with the vacuum energy density which 
does not dilute with the expansion pvac = A = const. ^ the energy 
density of CDM and baryons dilutes like oc a“^, while the density 
of radiation decays very quickly oc a~‘^. 

^ A de Sitter cosmological phase is a stage of the cosmological 
expansion whose dynamics is described by the line element ds^ = 
—dt^ where Hq is the Hubble constant. Since the 

deceleration parameter q = —1 — H/H^ = — 1 is a negative 
quantity, the de Sitter expansion is inflationary. 
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orbit joining these critical points then represents a feasi¬ 
ble image of the featured cosmic evolution in the phase 
space. 

While the theory of the dynamical systems has brought 
several surprises in the recent past in connection with 
the theory of chaos UM and fractals Q, the application 
of the tools of the dynamical systems in cosmological 
settings is less known in spite of the amount of published 
scientific pap ers on this subject (just to cite few of them 
see and the related references therein). 

This paper aims at an introduction to the application 
of the simplest tools of the theory of the dynamical sys¬ 
tems - those which can be explained with only previous 
knowledge of the fundamentals of linear algebra and of 
ordinary differential equations - within the context of 
the so called “concordance” or ACDM (section |TTO and 
scalar field models (section IIVI) of dark energy d, [n^ 
Emil- The dynamical systems tools play a central 
role in the understanding of the asymptotic structure of 
these models. I. e., these tools are helpful in looking 
for an answer to questions like: where does our Universe 
come from? or, what would be its fate? An specific ex¬ 
ample: the cosh potential, is explored in the section El 
in order to show the power of these tools in the search 
for generic dynamical behavior. 

We want to stress that there are very good introduc¬ 
tory and review papers on the application of the dynami- 
cal systems in cosmology - see, for instance, (a. 1 10l - ll2 . 1 16l| 
- so, what is the point of writing another introductory pa¬ 
per on the issue? In this regard we want to point out that, 
in the present paper, we aim not only at providing just 
another introductory exposition of the essentials of the 
subject but, also, at filling several gaps regarding specific 
topics usually not covered in similar publications exist¬ 
ing in the bibliography. Here we pay special attention, in 
particular, to the necessary (yet usually forgotten) defini¬ 
tion of the physically meaningful phase space, i. e., that 
region of the phase space which contains physically mean¬ 
ingful cosmological solutions. We also include the deriva¬ 
tion of relevant formulas in order to facilitate the way for 
beginners to concrete computations. We shall explain, 
in particular, a method developed in to consider ar¬ 
bitrary potentials within the dynamical systems study 
of scalar field cosmological models. This subject is not 
usually covered in the existing introductory bibliography. 
At the end of the paper, for completeness, we include an 
appendix section Em where an elementary introduction 
to the theory of the dynamical systems is given and use¬ 
ful comments on the interplay between the cosmological 
field equations and the equivalent phase space are also 
provided. Only knowledge of the fundamentals of linear 
algebra and of ordinary differential equations is required 


to understand the material in this section. The exposi¬ 
tion in the appendix is enough to understand how the 
computations in sections IIIII and IIVI are done. We rec¬ 
ommend those readers which are not familiar with the 
dynamical systems theory to start with this section. 

Our goal is to keep the discussion as general as pos¬ 
sible while presenting the exposition in as much as pos¬ 
sible pedagogical way. The paper contains many foot¬ 
notes with comments and dehnitions which complement 
the main text. Since the subject is exposed with some 
degree of technical details and a basic knowledge of cos¬ 
mology is assumed, this paper is suitable for teachers, 
undergraduate and postgraduate students from physics 
and mathematics disciplines. 


II. ASYMPTOTIC COSMOLOGICAL 
DYNAMICS: A MODEL-INDEPENDENT 
ANALYSIS 

In this section the asymptotic structure in the phase 
space of a general DE model will be discussed in detail. 
For simplicity and compactness of the exposition here 
we shall assume only a two-component cosmological fluid 
composed of cold dark matter - labeled here by the index 
“m,” i. e., {pm,Pm) {px,Px) - and of dark energy (the 
X-component). 

In order to derive an autonomous system of ordinary 
differential equations (ASODE) out of ([T]) it is useful to 
introduce the so called dimensionless energy density pa¬ 
rameters of matter and of the DE 


n 


m 


Pm 




Px 


(4) 


respectively, which are always non-negative quantities. 

In terms of these parameters the Friedmann equation - 
first equation in o - can be written as the following 
constraint: 


— 1) (5) 

which entails that none of the non-negative dimensionless 
energy density components alone may exceed unity: 0 < 
< 1, 0 < Hx < 1. 

Given the above constraints, if one thinks of the di¬ 
mensionless density parameters as variables of some state 
space, only one of them is linearly independent, say Hx. 
Let us write a dynamical equation for Hx, for which 
purpose we rewrite the Raychaudhuri equation - second 
equation in (ED - and the conservation equation for the 
X-component, in terms of 


For a very concise introductory knowledge of the theory of the 
dynamical systems and the related concepts like critical points 
(attractors and saddle points), heteroclinic orbits, etc, see the 
appendix [ym 


In agreement with conventional non-negativity of energy we con¬ 
sider non-negative energy densities exclusively. 
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2^2 — 3'Ym^m 

= -37m + 3(7^ - 7x)^x, 
^ = -SjxHnx. 

Then, if we substitute (O into 


rix 


PX 

3H^ 


— 2 — 


and consider the constraint ([S]), we obtain 


( 6 ) 


= 3(7m — 7x)f^x(l — f^x), (7) 

where the tilde denotes derivative with respect to the new 
variable r = Ina (12 = 7712', etc.). 

In general, for varying jx i another ordinary differential 
equation (ODE) for 'jx is needed in order to have a closed 
system of differential equations (the barotropic index of 
matter 7 ^ is usually set to a constant). The problem 
is that, unlike the ODE ([7]) which is model-independent, 
the ODE for 7 ^ = (px+Px)/Px requires of certain spec¬ 
ifications which depend on the chosen model (see below). 

Anyway, certain important results may be extracted 
from © under certain assumptions. For instance, if we 
assume a constant 7 x, the mentioned ODE is enough 
to uncover the asymptotic structure of the model © in 
the phase space. In this case the phase space is the I- 
dimensional segment 12x S [0,1]. The critical points are 


12x(I — 12x) = 0 12x = I, 12x = 0, 

where we are assuming that 7 x ^ 7 m- Tbe first fixed 
point 12 x = 1 is correlated with dark energy domination, 
while the second one 12x = 0 (12^ = I) is associated with 
CDM dominated expansion. 

If linearize Eq. © around the critical points: 1 — —>■ 
1, and 0 -I- do —>■ 0 (di and do are small perturbations), 
one obtains 


Besides, if 7 m > 7 x, the matter-dominated solution 
12m = 1 (l^X = 0): a{t) OC t2/37m^ is the past attractor 
since the initially small perturbation do(T) oc 
increases exponentially with r, thus taking the system 
away from the condition 12m = 1- Otherwise, if 7 m < 7 x, 
the X-dominated solution is the past attractor while the 
matter-dominated solution is the future attractor. How¬ 
ever, this last situation is not consistent with the known 
cosmic history of our Universe. 

The interesting thing about the above result is that 
it is model-independent, i. e., no matter which model 
to adopt for the DE (X-component), given that 7 x is a 
constant, the above mentioned critical points - together 
with their stability properties - are obtained. Besides, 
since a non-constant 7 x obeys an autonomous ODE of 
the general form 7 ^ = f{'yx,^x), at any equilibrium 
point where, necessarily 

12 x = 0 , 7 x = 0 ^ 7 x = 7 x = const. 

provided that jx Jm, either 12 x = 1 , or 12 x = 0 
(12m = !)■ Hence, independent of the model adopted to 
account for the X-component, matter-dominated and X- 
fluid dominated solutions are always equilibrium points 
in the phase space of the two-fluids cosmological model 
depicted by equations ©. This is a generic result. 


III. ACDM MODEL 

Let us now to show with the help of a specific cos¬ 
mological model of dark energy, how it is possible that 
a 2-dimensional system of ODE may come out of the 
3-dimensional system of second-order field equations o 
where, we recall, one of these equations is redundant. 
For this purpose we shall focus in the so called A-cold 
dark matter (ACDM) model whose phase space is a 
subspace of the phase plane.Here - only for the ex¬ 
tent of this section - in addition to the CDM and to the 
DE (the cosmological constant) we shall consider a radi¬ 
ation fluid with energy density pr and pressure Pr = Pr/S 
(pfr = 4/3). The conservation equation for radiation 
reads 


5[ = -i{im - lx)5i => ^ i ( t ) = 

<5^ = 3 ( 7 ™ - 7x)<Jo ^ 5 o(t) = 4 (8) 

As seen the DE-dominated solution: 

fix = I => 377^ = Px oc => a(t) oc 

is a stable (isolated) equilibrium point whenever 7 m > 
7 x since, in this case, the perturbation 

<5i(t) oc 

decreases exponentially. Hence, given that 7 m > 7 x, the 
dark energy dominated solution fix = I, is the future at¬ 
tractor solution which describes the fate of the Universe. 


Pr + AHpr = 0 => Pr- oc a 

Meanwhile, for the cosmological constant term - as for 
any vacuum fluid - we have that px = — Pa (7a = 0 ), 
which means that pA = 0 , so that this term does not 
evolve during the course of the cosmic expansion. 

In what follows, since we adopted the units system with 
SttG = c = I, we write px = A. Besides, since we deal 
here with CDM, which is modeled by a pressureless dust. 


For a dynamical systems study of the FRW cosmological model 
with spatial curvature fc 7 ^ 0 and with a cosmological constant 
see 0 . 
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we set 7m = 1- After these assumptions, and identifying 
px = Pa = A in equations ([T]), the resulting cosmological 
equations are 


where we have to take in mind the Raychaudhuri equa¬ 
tion and the conservation equation for radiation in (l9|) 
written in terms of the variables of the phase space 


3-ff^ — Pi- + Pm + A, H 


1 

2 



Pm “t“ — 0, Pr -\- AHPy- — 0. 


(9) 


Let us at this point to make a step aside to explain 
how it is that the cosmological constant can explain the 
present speed-up of the cosmic expansion. If we con¬ 
veniently combine the first two equations above, taking 
into account that straightforward integration of the last 
equations in (|9]) yield pm = GCmla^ and pr = SCrja^ 
respectively {Cm and Cr are constants), we get: 


^ = H + H^ 
a 



^ A 
a3 3' 


We see that acceleration of the expansion d > 0 may 
occur in this model thanks to the third term above 
(the DE, i. e., the cosmological constant). Besides, 
the accelerated pace of the expansion is a recent phe¬ 
nomenon: during the course of the cosmic evolution the 
initially dominating radiation component dilutes very 
quickly oc a~'^ until the matter component starts dom¬ 
inating. As the cosmic expansion further proceeds the 
matter component also dilutes oc a~^ until very recently 
when the cosmological constant started to dominate to 
yield to positive d > 0. 

Following the same procedure as in the former section 
we write the Friedmann constraint 


H = ——{3 + x-3y), 

and then we have to go to derivatives with respect to the 
variable r = Ina(t) (^' —We obtain: 


x'=-x{l-x + 3y), y'= {3 + x-3y)y. (13) 

It is remarkable the simplicity of the system of two 
ordinary differential equations m as compared with 
the system of three second-order cosmological equations 
®. The critical points of (fT51) in 4 'a are easily found 
by solving the following system of algebraic equations: 
a: (1 — a: -I- 3y) = 0, (3 -I- a: — 3y) y = 0. 

In the present case, in terms of the variables x, y, the 
deceleration parameter q = —I — H/H'^ (equation ([3])) 
can be written as: 


q= {l + x-3y)/2, 

which means that those fixed points which fall in that 
region of the phase space laying above the line y = 
1/3 -I- x/3, correspond to cosmological solutions where 
the expansion is accelerating. 


A. Critical points 


Hj. -i- Urn + — 1) (10) 

where fir = pr/3H^ and Ha = A/3H^. In what follows 
we choose the following variables of the phase space: 

X = Hri y — HAi (11) 

where, for sake of simplicity of writing, we adopt x and y 
in place of the dimensionless energy densities of radiation 
fir and of the cosmological constant flA respectively. But 
we warn the reader that the same symbols x and y will 
be used in the next sections to mean different variables 
of the phase space. 

We have that fl^ = ^ — x — y, and, since 0 < a; < 1, 
0 < y < 1, and 0 < flm < 1, the physically relevant phase 
space is defined as the following 2D triangular region (see 
the top panel of FIG. [J): 

^'a = {ix,y) : a:-|-y<l, 0<a;<I, 0<y<l}. (12) 

As before, in order to derive the autonomous ODE-s, 
we first take the derivative of the variables x and y with 
respect to the cosmic time t: 

. _ Pr n, Pr H . ^ ^ A H 

'^~3m 3H^H^y~ 3m 


1. Radiation-dominated critical point 

The critical point Pr : (1 j0)j fir = 1 => 3iL^ = pr 
a{t) oc v^, corresponds to the radiation-dominated 
cosmic phase. This solution depicts decelerated expan¬ 
sion since 9=1. The linearization matrix for the system 
of ODE (fT31) - for details see the appendix - is 

— l-\-2x — 3y —3x \ 

y 3 X — 6yJ ' 

The roots of the algebraic equation 

det I J(Pr) - API = det ^ = 0, 

are the eigenvalues Ai = 1 and A 2 = 4. Since both eigen¬ 
values are positive reals, then Pr is a source point (past 
attractor). This means that in the ACDM model de¬ 
scribed by the cosmological equations ® the radiation- 
dominated solution a{t) (X ^/t is a. privileged solution. 

From the inspection of the phase portrait of (fT3ll - top 
panel of FIG. [1] - it is apparent that any viable pattern 
of cosmological evolution should start in a state where 
the matter content of the Universe is dominated by radi¬ 
ation. Of course this is a drawback of our classical model 
which can not explain the very early stages of the cosmic 
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FIG. 1: Phase portrait of the ASODE (I13II corresponding to 
the ACDM model (top panel). The orbits of (1131) emerge 
from the past attractor Pr : (1,0) - the radiation domination 
state - and end up in the future attractor Pds '■ (0,1) (de 
Sitter phase). For a non-empty set of initial conditions the 
corresponding orbits approach close enough to the matter- 
dominated saddle point at the origin Pm : (0,0). In the bot¬ 
tom panel the phase portrait of the ASODE m - showing 
two orbits generated by two different sets of initial conditions 
- is shown. Here the free parameter jj, has been arbitrarily 
set fj, — 10. 


evolution where quantum effects of gravity play a role. 
This includes the early inflation. In a model that would 
account for such early period of the cosmic dynamics, 
inflation should be the past attractor in the equivalent 
phase space. 


2. Matter-dominated critical point 

The critical point Pm ■ (0,0), flm = 1 => = pm 

=> a{t) oc corresponds to the matter-dominated 
solution, which is associated with decelerated expansion 
{q= 1/2). Following the same procedure above - see the 
appendix I VIII - one finds the following eigenvalues of the 
linearization matrix J{Pm) evaluated at the hyperbolic 
equilibrium point Pm'- Ai = 1, A 2 = —3. Hence this is 
a saddle critical point. As seen from the inspection of 
the phase portrait in the top panel of FIG. [U only for 


conveniently chosen initial conditions the given orbits in 
'Fa approach close enough to Pm- Since this is a unstable 
(metastable) critical point, it can be associated with a 
transient stage of the cosmic expansion only. This is good 
since a stage dominated by the dark matter can be only a 
transient state, lasting for enough time as to account for 
the observed amount of cosmic structure. A drawback is 
that we need to fine tune the initial conditions in order 
for feasible orbits in the phase space to get close enough 
to Pm- 


3- de Sitter phase 

The critical point Pds ■ (0,1), Ha = 1 3H^ = A 

=> a{t) oc corresponds to a stage of inflationary 

(q = — 1) de Sitter expansion. The eigenvalues of the ma¬ 
trix J{Pds) are Ai = —3 and A 2 = —4. This means that 
Pds is a future attractor (see the classification of isolated 
equilibrium points in the appendix IVIIl) . This entails that 
independent on the initial conditions chosen H(!, H((, the 
orbits in IFa are always attracted towards the de Sit¬ 
ter state which explains the actual accelerated pace of 
the cosmic expansion in perfect fit with the observational 
data. This is why, in spite of the serious drawback in con¬ 
nection with the cosmological constant problem 111 , 123 , 
the very simple ACDM model represents such a success¬ 
ful description of the present cosmological paradigm and 
it is, therefore, coined the “concordance model”. Any 
other DE model, regardless of its nature, ought to be 
compared to the ACDM predictions as a first viability 
test. 

Additional refinement of the above model is achieved 
if we add the energy density of baryons. 


IV. SCALAR FIELD MODELS OF DARK 
ENERGY: THE DYNAMICAL SYSTEMS 
PERSPECTIVE 

The cosmological constant problem can be split into 
two questions [21]: (i) why the vacuum energy p-^ac = A 
is not very much larger? - old cosmological constant 
problem - and (ii) why it is of the same order of magni¬ 
tude as the present mass density of the universe? which is 
acknowledged as the new cosmological constant problem. 
In order to avoid the old cosmological constant problem, 
which is exclusive of the ACDM model, scalar field mod¬ 
els of dark energy are invoked.In this last case an ef¬ 
fective “dynamical” cosmological constant is described by 
the scalar field’s X self-interaction potential V{X). This 
may be arranged in such a way that at early times the 
vacuum energy pvac.o = V{Xq) {X = 0) is large enough 


The new cosmological problem is inherent also in several scalar 
field models of dark energy. 
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as to produce the desired amount of inflation, while at 
late times Pvacj = V{Xf) is of the same order of the 
CDM energy density pm- 

In this section, we shall explore the asymptotic struc¬ 
ture of general scalar field models of DE [ 3 , llll - fl^ [H- 
[12) 113) 112] • These represent a viable alternative to the 
ACDM model explored above. Among them “exponen¬ 
tial quintessence” V{X) = Mexp{—pX) 0111 is one 
of the most popular scalar field models of dark energy. 

The specification of a scalar field model for the DE 
means that the energy density and the parametric pres¬ 
sure in the cosmological field equations o are given by 


One of these constraints is given by 

0 < fim < 1 ^ 0 < + y'^ < 1. 

Besides, since 0 < fljc < 1, then \x\ < 1. Additionally we 
shall be interested in cosmic expansion exclusively H > 0, 
so that y > 0- Hence, for instance, for potentials of one of 
the following kinds (see below): V = Vq, V = 
the physical phase plane is defined as the upper semi¬ 
disk, 

:= {(a;, y) : 0 < -f < 1, |a;| < 1, y > 0}. (19) 


px=X‘^l2 + V{X),px=X^l2-V{X), (14) 

respectively. In these equations V = V{X) is the self¬ 
interacting potential of the scalar field X. Besides, for 
the EOS parameter ujx = 7x — 1, one has 


However, as we shall see, in general - arbitrary potentials 
- the phase space is of dimension higher than 2 [^ . 


B. Autonomous system of ODE-s 


px A2_^21/' 


(15) 


After the above choice, the conservation equation © can 
be written in the form of the following Klein-Gordon 
equation: 


X + -iHX = -V,x- (16) 

It happens that deriving an autonomous ODE for the 
variable ojx - see section |H]- can be a very difficult task. 
Hence it could be better to choose a different set of phase 
space variables which do the job. Here, instead of the 
phase space variables fix and ujx we shall choose the 
following variables 0: 


In order to derive the autonomous ordinary equations 
for the phase space variables x, y one proceeds in a similar 
fashion than in section im First we writte the Raychaud- 
huri equation in © 


^ — 2 (TiriPm Px Vx) ^ 

H 

-2-^ ~ , 
in terms of the new variables: 


- 2^ = 37m (1 - - y^) + 6a:^, (20) 

where we have taken into account the Friedmann con¬ 
straint ®. Then, given the definition dm, let us find 


= ^ 


(17) 


so that the dimensionless energy density of the scalar 
field X is given by Hx = x'^ + y"^, and the Friedmann 
constraint © can be written as 


X X H 

Vgh Vgh h ’ 


( 21 ) 


or if one substitutes m and (1201) back into (1211) . then 


Hm = 1 - - y^. 


(18) 


x' = —3a:(l — x^) + 



A. Determination of the physical phase space 

The first step towards a complete study of the asymp¬ 
totic structure of a given cosmological model is the rig¬ 
orous determination of the phase space where to look for 
the relevant equilibrium points. In the present model 
there are several constraints on the physical parameters 
which help us to define the physically meaningful phase 
space. 


where we replaced derivatives with respect to the cosmic 
time t by derivatives with respect to the variable r = 
Ina(t). Applying the same procedure with the variable y 
one obtains 

, flVxX h\ 

^ V H 

where we have taken into account that V = V,xX, hence 









Correspondence f{z) V{X) 


Function f{z) 

z 

Self-interaction potential V{X) 

Reference 

fiz) = 0 

±y 

V = Vo 

mi 

f{z) = -al3 -1- (a -1- I3)z - z^ 

a+P 

l+e(/3-c)X 

U = Uo (e“-^ -b e'^^) 


fiz) = (pV" - z^)/p 

py cotanh(/iA) 

V = Vo sinhP(pA) 


fiz) = - z^)/p 

pptanh(pA) 

V = UocoshP(pA) 


f{z) = (pV" - z^)/2p 

sinh(nX) 
coshf/LX) —1 

V = Vo [cosh(pA) - 1]^ 

[32] 


TABLE I: The functions f{z) in Eq. (I27II for different self-interaction potentials. 


y' = y 



y^) + 3a;^ 



(23) 


Unless V = Vo ^ Vx/V = 0, or U = Vb => 
Vx/V = ±/r = const., the ASODE (|2^ . (|2^ is not a 
closed system of equations since one equation involving 
the derivative of Vx /U with respect to r is lacking. 

One example where the ASODE (j^ . (1^ is indeed 
a closed system of ODE-s is given by the so called “ex¬ 
ponential quintessence.” This case has been investigated 
in all detail in m so that we submit the reader to that 
reference to look for a very interesting example of the ap¬ 
plication of the dynamical systems tools in cosmology. As 
a matter of fact, it results a very useful exercise for those 
who want to learn how to apply the dynamical systems 
tools in cosmological settings to reproduce the results of 
the dynamical systems study in the seminal work [HI- 

In order to be able to consider self-interaction poten¬ 
tials beyond the constant and the exponential potentials, 
in addition to the variables x, y, one needs to adopt a 
new variable [3, [3, 


z = VxlV, (24) 

so that z = 0 corresponds to the constant potential, while 
z = ±/r is for the exponential potential.^'* By taking the 
derivative of the new variable z with respect to r one 
obtains 


z' = V6xf{z), (25) 

where we have defined (as before take into account that 
V = VxX and Vx = Vxx^, etc.) 


For non-exponential potentials, besides the increase in the di¬ 
mensionality of the phase space, there is another problem. Typ¬ 
ically the phase space becomes unbounded, so that two different 
sets of variables are required to cover the entire phase space 
|20ll . In the present paper, however, the chosen example of the 
cosh-like potential does not lead to unbounded phase space. In 
consequence a single set of phase space variables is enough to 
describe the whole asymptotic (and intermediate) dynamics. 


f{z) = z^T-l],r = VVxx/V}„ (26) 

and the main assumption has been that the above E is a 
function of the variable z: F = r(z). Notice that, since 

UUxx _ V.xxjV _ 1 Uxx 

y.x (Vx/V)^ z^ U ’ 

the left-hand equation in (1261) can be rewritten also as 


f{z) = ^-z\ (27) 

In the event that F can not be explicitly written as a 
function of z, then, an additional ODE: F' = ..., is to be 
considered. However, this case is by far more complex 
and does not frequently arise. 

Before going further we want to point out that, as a 
matter of fact, the function F and the variable 2 were 
first identified in [s^. Besides, the dynamical system 
(1^ , (1^ , (1^ was explored in M several years before 
it was studied in the reference [23|. However, it was in 
[3 | where the possibility that for several specific self¬ 
interaction potentials F can be explicitly written as a 
function of z, was explored for the first time. In the ref¬ 
erence [T3 |- although the correct dynamical system was 
identified, the authors were not interested in specific self¬ 
interaction potentials. The cost of the achieved general¬ 
ity of the analysis was that the authors had to rely on 
the obscure notion of “instantaneous critical points”. 


V. AN EXAMPLE: THE COSH-LIKE 
POTENTIAL 

In order to show with a concrete example how the func¬ 
tion f{z) can be obtained for an specific potential, let us 
choose the cosh potential [3 HI] 


U = Uocosh(/rA). (28) 

This potential - as well as a small variation of it in the 
last row of TAB. |T]- has a very interesting behavior near 
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of the minimum at X = 0. Actually, in the neighborhood 
of the minimum, the potential (I28p approaches to 

y(X) « tb + m2 = 

The quadratic term is responsible for oscillations of the 
scalar field around the minimum which play the role of 
cold dark matter, hence, at late times when the system 
is going to stabilize in the minimum of V{X) the cosh 
potential makes the scalar field model of DE to be in¬ 
distinguishable from the standard ACDM model (see the 
properties of the equivalent model of Ref. El)- 

Taking derivatives of (1281) with respect to the scalar 
field X we have that Vx = fJ,Vosmh(pX), and Vxx = 
Hence, 


z = y.xjy = ^tanh(^X), f{z) = fj,'^ — (29) 

Notice that the z-s which solve f{z) = 0, i. e., z = 
±fj,, which correspond to the exponential potentials, are 
critical points of (HH). In TAB.U the functions /(z) for 
several well-known potentials are displayed. 

Let us to collect all of the already found autonomous 
ODE-s which correspond to a scalar field model of DE 
with arbitrary self-interaction potential: 


x' = —3x(l — a;2) -|- 



z' = V6xf{z), (30) 

where it is evident that the equations for x and y are 
independent of the self-interaction potential, and that 
details of the given model are encoded in the function 
/(z) - Eq. (imi ~ which depends on the concrete form of 
the potential (see TAB.|I|. 


A. Critical points 

If we apply the simplest tools of the dynamical systems 
theory to the the above example with the cosh poten¬ 
tial the corresponding asymptotic dynamics in the phase 
space is revealed. Given that, at the minimum of the cosh 
potential, the model behaves as ACDM, here we shall as¬ 
sume that the matter fluid is composed mainly of baryons 
so that it behaves like pressureless dust {'jm = !)• 

The following closed system of autonomous ODE-s is 
obtained (it is just Eq. (1501) with the appropriate substi¬ 
tutions): 



Since z = /rtanh(/xX) => —y < z < fi, the phase space 
where to look for critical points of (lOTll is the bounded 
semi-cylinder (see Eq. (HU)): 

^'cosh := {ix,y,z) : 0 < x"^ + y'^ < I, 

\x\ < 1, y > 0, |z| < y}. (32) 

Nine critical points Pi : {xi,yi,Zi) of the autonomous 
system of ODE (1511) and one critical manifold are found 
in ikcosh- These together with their main properties are 
shown in TAB. im Eor sake of conciseness, here we shall 
concentrate only in a few of them which have some par¬ 
ticular interest. 

The de Sitter critical point : (0,1,0), for which 
Ox = 1 (0„ = 0), is associated with accelerating expan¬ 
sion since the deceleration parameter 

iL 1 3 , 9 o, 

9 = -1-;^ = 2 + 2^"^ 

is a negative quantity. Since z = 0 => E = Vb, and 
since a; = 0 ^ A = 0, and 

y = l ^ H = ^ a(t) oc e^*, 

this point corresponds to a de Sitter (inflationary) phase 
of the cosmic expansion. We underline the fact that to a 
given point in the phase space - in this case Pds G ^cosh “ 
it corresponds a specific cosmological dynamics (de Sitter 
expansion in the present case). 

Next we have the scalar field A-dominated equilibrium 
points 

The first thing one has to care about is the existence of 
the critical point. By existence we mean that the given 
point actually belongs in the phase space ('Lcosh in the 
present case). Hence, since |x| < 1 => < 6- This 

is precisely the condition for existence of the points P^: 
y? < 6. This warrants, besides, that y = y^l — be 
a real number. If one substitutes (ED back into ED on 
gets 

2 2 
X — y 
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Critical Points 

X 

y 

z 

Existence 

0 ,x 


Ai 

A2 

As 

Z 

0 

0 

z 

Always 

0 

1 

3/2 

-3/2 

0 

PdS 

0 

1 

0 

?? 

1 

0 

_3 , 1 Jl _ TZ . 

2 ^ 2 V 3 

3 3/1 4^2 

2 2 V 3 

-3 


±1 

0 

d 

5? 

1 

0 

3±^M 

=F2%/67i 

3 


±1 

0 


?? 

1 

0 

3 T ■\J^P 

±2V6p 

3 

p± 

X jm 

V ^ ^ 

/¥l 

y 2 fj , 

TP 

I 

3 

1-4 

3,3 / 24 V 

-4 + 

3 Z j 24 . 7 

4 4V7^ 

6 

p± 

± JL . 


TP 

p ^ <6 

1 

0 

- iTp ^ 

-3+4 

2 p ^ 


TABLE II: Critical points of the autonomous system of ODE (1311) and their properties. 


SO that, for the present case the EOS parameter of the X- 
field ujx = 1 7x = 1*^/3. Then the continuity 

equation o can be readily integrated at to obtain: 
px = M = M a~^ , where M is an integration 

constant. After this the Friedmann constraint 

Ox = 1 ^ = px ^ 

a V 3 

can be integrated also to obtain the following cosmolog¬ 
ical dynamics a(t) oc , which is to be associated 
with the critical points P^. We want to point out that 
whenever 2 < < Q {q > 0) the present scalar field 

dominated solution represents decelerated expansion. 

In the case of the matter-dominated critical manifold^® 

M = {{0,0,z):z€R}, 0 „ = 1 , <7 = 1 / 2 , 


one needs to linearize the system of ODE - in this case 
- which means, in the end, that we have to find the 
eigenvalues of the linearization or Jacobian matrix J (see 
Eq. (1551) of the appendix) 


J = 


^dx'jdx dx'/dy dx'/dz\ 
dy'/dx dy'/dy dy'/dz 
\dz'jdx dz'jdy dz'jdzj 


Next J is to be evaluated at each non-hyperbolic critical 
point, for instance, J{Pds), J{Px), etc. The eigenvalues 
of the resulting numerical matrices is what we use to 
judge about the stability of the given points. 

Take, for instance, the de Sitter point Pds '■ (0,1,0). 
In this case one has to find the roots of the following 
algebraic equation for the unknown A: 


instead of an isolated equilibrium point this is actually 
a 1 -dimensional manifold which is extended along the z- 
direction. This means in turn that the matter-dominated 
solution: = 1 


det 


/-3 - A 0 

0 -3- A 0 

\-y/Qp 0 -A 


= 0 . 


3iL^ = Pm = Ma ^ => a{t) oc 


The resulting eigenvalues are 


may coexist with the scalar field domination no matter 
whether V oc cosh( 7 iX) —intermediate X-s, V = Vq ^ 
AT = 0, or E oc —>■ very large X-s (|X| —>• oo). 

This entails that it can be found an heteroclinic orbit in 
the phase portrait - an orbit connecting two or more dif¬ 
ferent equilibrium points - joining the matter-dominated 
equilibrium state with a scalar field domination critical 
point, be it originated either by a cosmological constant, 
or by exponential or cosh-like potentials. 

B. Stability 

In order to judge about the stability of given hyperbolic 
equilibrium points (as explained in the appendix Em 


Here by critical manifold we understand a curve in the phase 
space all of whose points are critical points. 


Xp2=-i±yi-^p^, a3=-3. 

For < 3/4 the critical point is an isolated focus, while 
for p'^ > 3/4, since the eigenvalues Ai and A 2 are com¬ 
plex numbers with negative real parts, the point Pds is 
a stable spiral. Hence the de Sitter solution is always a 
future attractor (see the bottom panel of FIG. [!])• This 
is consistent with the fact that the cosh potential is a 
minimum at AT = 0 => V = Vq. 

The eigenvalues of the linearization matrix J{P^) are 

Ai = —3 -l- p^/2, X 2 = —3 -l- p^, A 3 = p^. 

Hence, since as required by the existence p'^ < 6, P^ are 
always saddle critical points (at least one of the eigenval¬ 
ues is of a different sign). 

If one computes the eigenvalues of J{M) one finds 
Ai ^2 = ±3/2, A 3 = 0. The vanishing eigenvalue is due to 
the fact that the 1 -dimensional manifold extends along 



























11 


the 2 -direction.^® The differing signs of the eigenvalues 
Ai and A 2 means that each point in Ad is a saddle. It 
might be a nice task for the interested readers to judge 
about the stability of the remaining critical points of m 
(see TAB. HU). 


of Mexico. Last but not least the authors are indebted 
to the anonymous referees for the constructive criticism. 

VII. APPENDIX: DYNAMICAL SYSTEMS IN 
COSMOLOGY 


VI. CONCLUSION 

In this paper we have shown in as much as possible 
pedagogical way, how the application of the basic tools 
of the dynamical systems theory can help us to deeper 
understand the dynamics of cosmological models of dark 
energy. Our main aim has been to make clear how much 
useful information about the asymptotic cosmological dy¬ 
namics - the one which decides the origin and the fate of 
our Universe - can be extracted by means of these simple 
tools. We have concentrated our exposition in scalar field 
cosmological models of dark energy and the analysis was 
kept as general as possible. Special attention has been 
paid to the derivation of relevant equations and formu¬ 
las, so as to facilitate the way for beginners to concrete 
computations. 

It has been demonstrated that certain generic asymp¬ 
totic behavior arises which is quite independent on the 
concrete model of DE: matter dominance transient phase 
and dark energy dominance at late times. Just for illus¬ 
tration we have chosen a concrete potential - the cosh 
potential - which at late times approaches to the ACDM 
model. The approach undertaken here can be easily ap¬ 
plied to other self-interaction potentials (see TAB. HI. 

In the present paper, for completeness, we have in¬ 
cluded a very concise and simple exposition of the fun¬ 
damentals of the theory of the dynamical systems - see 
the appendix IVIII - which may be considered as material 
for an introductory course on this subject. 
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In general a critical point for which at least one of the real parts 
of the eigenvalues of the Jacobian matrix is vanishing, is called 
as a non-hyperbolic critical point. In such case the Hartman- 
Grobman theorem [jSl can not be applied. As a consequence 
one has to go beyond the linear stability theory. The centre 
manifold theory is useful in this case. For the application of 
the centre manifold theory in a case of cosmological interest see 

mm- 


In classical mechanics the state of a given physical sys¬ 
tem composed, say, of N particles, is completely speci¬ 
fied by the knowledge of the N generalized coordinates Qi 
{i = 1,..., N) and the N conjugated momenta pi, which 
satisfy the Hamilton’s canonical equations of motion 

= dU/dpi, Pi = -dH/dq^, 

where % = 'H{qi, ...,qM,Pi, ■■■,Pn jt) is the Hamilto¬ 
nian of the physical system. Hence, a given state of 
this system corresponds to a point in the 2N- dimen¬ 
sional space of states spanned by the 2N coordinates 
qi, ...,qN,Pi, ■■■,PN, also known as the “phase space.” 

Like it is for any other physical system, the possible 
states of a cosmological system may be also correlated 
with the points in an equivalent state space. However, 
unlike in the classical mechanics case, the choice of the 
variables - generalized coordinates and their conjugate 
momenta - for a cosmological model is not a trivial issue. 
In this case a certain degree of uncertainty in the choice of 
an appropriate set of variables of the phase space arises. 
There are, however, certain not written rules one fol¬ 
lows when choosing appropriate variables of the phase 
space: (i) the variables should be dimensionless, and (ii) 
whenever possible, these should be bounded. The latter 
requirement is necessary to have a bounded phase space 
where all of the existing equilibrium points are “visible,” 
i. e., none of then goes to infinity. 

The interplay between a cosmological model and the 
corresponding phase space is possible due to an existing 
one-to-one correspondence between exact solutions of the 
cosmological field equations (ED and points in the phase 
space spanned by given variables x, y}'^ 

X = x{H, p^, px), y = y{H, p^, px). 

When we replace the original field variables H, p^, and 
Px, by the phase space variables x, y, we have to keep in 
mind that, at the same time, we trade the original set of 
non-linear second order differential equations in respect 
to the cosmological time t (equations ED)) by a set of 
first order ordinary differential equations (ODE): 


x' = f{x, y), y' = g{x, y), (34) 

where the tilde denotes derivative with respect to the 
dimensionless “time” parameter, r = j da ja => dr = 


Here, for simplicity of the exposition, we consider a 2D phase 
space. 
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din a = Hdt. The introduction of r instead of the cosmic 
time t is dictated by simplicity and compactness of writ¬ 
ing. Besides, this warrants that the first order differential 
equations (IMl) are a closed system of equations. 

The most important feature of the system of ODE (IMl) 
is that the functions f{x,y) and g{x,y) do not depend 
explicitly on the parameter r. This is why (|M| is called as 
an autonomous system of ODE. The image of the integral 
curves of (IMl) in the phase space are called “orbits” of the 
system of autonomous ODE. The critical points of (p^ - 
also fixed or stationary points - Per : (a^ender) are those 
for which the “velocity” vector v = {x',y') vanishes (see 
below): 

v(Pcr) = {x'{Per), y'{Per)) = 0. 

In order to judge about the stability of given equi¬ 
librium points one needs to linearly expand (1341) in the 
neighborhood of each hyperbolic critical point: 

f{Xer + 5x, yer + Sy), g{Xer + Sx, J/er + Sy). 

Depending on whether the linear perturbations Sx = 
Sx{t), Sy = Sy{T) decay (grow) with time r or decay in 
one direction while grow in the other, the critical point 
can be a future (past) attractor or a saddle point. 

Knowledge of the critical - also equilibrium or fixed 
- points in the phase space corresponding to a given 
cosmological model is a very important information. In 
particular, the existence of attractors can be correlated 
with generic cosmological solutions that might really de¬ 
cide the fate and/or the origin of the cosmic evolution. 
The “meta-stable” saddle equilibrium points - station¬ 
ary points which are not local extrema - can be asso¬ 
ciated with (not less important) transient cosmological 
solutions. Hence, in the end we are trading the study 
of the cosmological dynamics depicted by P = H{t), 
Px = Px{t), px = Px{t), by the study of the properties of 
the equilibrium points of the equivalent autonomous sys¬ 
tem of ODE (IMl) in the phase space d' of the cosmological 
model, which, we recall, differs from the usual definition 
in classical mechanics. 


A. Remarks on phase space analysis 

In this subsection we provide a very simplified exposi¬ 
tion of the fundamentals of the dynamical systems the¬ 
ory which are useful in most cosmological applications 
@1 [m, [13 ■ Only knowledge of elementary linear algebra 
and of ordinary differential equations is required to un¬ 
derstand the material exposed here. Eor a more formal 
and complete introduction to this subject we recommend 
the well-known texts [1, [13 • 

The critical points of the system of ODE (l34l) Pi : 
{xi\yi),\. e., the roots of the system of algebraic equa¬ 
tions 

f{x,y) = 0, g{x,y) = 0, 


correspond to privileged or generic solutions of the orig¬ 
inal system of cosmological equations. To judge about 
their stability properties it is necessary first to linearize 
(1341) around the hyperbolic equilibrium points.^® To lin¬ 
earize around a given hyperbolic equilibrium point Pi 
amounts to consider small linear perturbations 

X ^ Xi + Sx{t), y ^ yi + Sy{T). 

These perturbations would obey the following system of 
coupled ODE (here we use matrix notation): 

<5x' = J{Pi) ■ <5x, <5x = i) ■ 

where J - the Jacobian (also linearization) matrix, is to 
be evaluated at Pi. Thanks to the Hartman-Grobman 
theorem [^, which basically states that the behavior 
of a dynamical system in the neighborhood of each hy¬ 
perbolic equilibrium point is qualitatively the same as 
the behavior of its linearization, we can safely replace 
the study of the dynamics of (l34l) by the corresponding 
study of its linearization (|55|) . 

We assume that J can be diagonalized, i. e., Jd = 
where M is the diagonalization matrix and 



is the diagonal matrix whose non vanishing components 
are the eigenvalues of the Jacobian matrix J: 

det\J-XU\=0,U=(^^^^. (36) 

After diagonalization the coupled system of ODE (1551) 
gets decoupled: 

Sx' = Jd ■ Sx, Sx = M~^Sx, (37) 

where we have to recall that to each equilibrium point it 
corresponds a different matrix Jd- The decoupled system 
of ODE ((571) is easily integrated: 

Sx{t) = Sx{0)e^^'', Sy{T) = Sy{0)e^^''. 

Since the diagonal perturbations Sx and Sy are linear 
combinations of the perturbations Sx, Sy: 

Sx = ciiSx + ci 2 Sy, Sy = C 2 iSx + C 22 Sy, 


Ph is an hyperbolic critical point if the real parts of all of the 
eigenvalues of the linearization matrix around are necessarily 
non-vanishing. In particular, an hyperbolic point can not be a 
center. In more technical words: an hyperbolic critical point is 
a fixed point that does not have any center manifolds. 
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FIG. 2: Phase portrait: Isolated equilibrium points in the phase space when the eigenvalues Ai and A 2 are real numbers (top 
panels). The velocity field, represented by the arrows field in the figures, shows where the flux of the ODE (1341) flows. Several 
orbits generated by different sets of initial conditions [2:(0) = xoi\y{Q) = yoi\ are shown in each case. In the left-hand top panel 
the origin P : (0; 0) is a unstable node or past attractor (A 2 > Ai > 0), while in the center top panel the origin is a saddle 
critical point (Ai < 0 < A 2 ), and in the right-hand top panel P : (0; 0) is a stable node or future attractor (Ai < A 2 < 0). In 
the bottom panels the equilibrium points for the case when the eigenvalues are complex numbers [Ai ,2 = iuj] are shown. In 
the left-hand bottom panel the origin (a;; y) = (0; 0) is a unstable spiral (i/ > 0), while, in the center panel it is a center {v = 0). 
The point (0; 0) is a stable spiral in the right-hand bottom panel {v < 0). 


where the constants Cij are the components of the matrix 
M~^, then 

Sx^t) = cii + C 12 

5 t/(T) = C21 -h C22 ( 38 ) 

where 

C22Sx{0) _ ci2(5y(0) 

Cll = -, C 12 =-, 

C22C11 — C12C21 C22C11 — C12C21 

C2i5x{ 0) _ cii5y{0) 

C 21 =-, C 22 =-. 

C22C11 - C12C21 C22C11 - C12C21 

As a matter of fact we do not need to compute the 
coefficients c^, instead the structure of the eigenvalues 
Ai is the only thing we need to judge about the stability 
of given (hyperbolic) equilibrium points of (l34l) . If the 
eigenvalues are complex numbers A± = u ± iuj the per¬ 
turbations (|Mll do oscillations with frequency uj. If the 
real part v is positive the oscillations are enhanced, while 
if < 0 the oscillations are damped. The case ^ = 0 is 
associated with a center - harmonic oscillations - and is 
not frequently encountered in cosmological applications. 
This latter kind of point in the phase space - eigenvalue 


with vanishing real part - is called as non-hyperbolic crit¬ 
ical point. 

B. Taxonomy of isolated equilibrium points 

Here we summarize the basic classification of isolated 
equilibrium points in the phase plane. For simplicity, we 
assume that Ai < A 2 , but nothing changes if one assumes 
that Ai > A 2 , or if one makes no assumption at all. In 
the figure!^ the illustrative phase portraits - the drawing 
of the trajectories of a dynamical system in the phase 
plane - are shown. 

1. The eigenvalues Ai, A 2 in Eq. (1551) are real num¬ 
bers. 

(a) 0 < Ai < A 2 - the critical point Pi is an unsta¬ 
ble node or, also, a source point (past attrac¬ 
tor). This is a unstable critical point which 
represents the origin of a non-empty set of or¬ 
bits in the phase plane (left-hand top panel of 
FIG. ID). 

(b) Ai < 0 < A 2 - Pi is a saddle point. This is 
a unstable equilibrium point but can be as¬ 
sociated with a marginally stable state since 















14 


the orbits of the system of autonomous ODE-s 
(IMl) spend some “time” in the neighborhood 
N{e,Pi) of the critical point until these leave 
N{e,Pi) to go elsewhere in the phase plane 
(center top panel of FIG. [2]). 

(c) Ai < A 2 < 0 (also Ai = A 2 < 0) - Pi is a stable 
node (future attractor). This kind of equilib¬ 
rium point is associated with an asymptoti¬ 
cally stable state. The future attractor is the 
end point of a non-empty set of phase space 
orbits generated by a wide range of initial con¬ 
ditions (right-hand top panel of FIG. [2]). 

2. The eigenvalues Ai, A 2 in Eq. (l38l) are complex 
numbers: Ai ,2 = v ±iui (co ^ 0). 

(a) > 0 - the equilibrium point Pi is an unstable 

spiral (past attractor). It is associated with 
unstable oscillations, i. e., with amplified os¬ 
cillations (left-hand bottom panel in FIG. [5]). 


(b) 1 / = 0 - Pi is a center (not frequent). It is as¬ 
sociated with free oscillations (center bottom 
panel in FIG. [2]). 

(c) v <0 - Pi is a stable spiral (future attractor). 
The system does damped oscillations until it 
settles down in the equilibrium state (right- 
hand bottom panel in FIG. [2]). 


The above classification of isolated critical points en¬ 
compasses the kinds of equilibrium points most fre¬ 
quently encountered in cosmological applications. Any 
curve in the phase space which connects critical points 
(two or more, perhaps all of them) is called as “hetero¬ 
clinic orbit.” A curve which joints a critical point with 
itself is called “homoclinic trajectory.” There are other 
useful concepts which are not explained here. The inter¬ 
ested reader is submitted to the bibliography BEIli- 
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